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Abstract 

A possibility of the electronic origin of the high-temperature superconductiv- 
ity in cuprates is probed with the quantum Monte Carlo method by revisiting 
the three-band Hubbard model comprising Cu3dx2_y2 and 02p„ orbitals. The 
d^2_y2 pairing correlation is found to turn into an increasing function of the 
repulsion Ud within the d orbitals or the d-p level off-set Ae, where the corre- 
lation grows with the system size. We have detected this in both the charge- 
transfer and Mott-Hubbard regimes upon entering the strong-correlation re- 
gion (Ud or Ae > bare band width). 
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The high-temperature superconductivity in cuprates harbors some of the most fascinat- 
ing aspects of strongly-correlated electron systems, but, despite a body of theoretical works 
on possible electronic mechanisms of superconductivity, we are still someway from a com- 
plete understanding of what happens in the realistic parameter range. Experimental and 
theoretical studies have indicated that the essence of the cuprates lies in the two-dimensional 
Cu02 plane, for which it is generally recognized that Emery's three-band Hubbard model 
im is the basic, starting model that describes both the copper 3d and oxygen 2px and 2py 
orbitals. 

The model captures the essential feature of the system with two key parameters: Ud 
(the on-site Coulomb repulsion between copper d holes) and Ae {Cu3d-02p level offset), 
where the energies are measured in units of the d-p hybridization, t^p. The presence of both 
parameters makes the physics richer. Specifically, the inequality Ae < Ud is usually used 
to identify the insulating host material as a charge-transfer insulator, as opposed to the 
Mott-Hubbard insulator with Ae > Ua^- Here we shall extend this terminology into the 
doped case. The three-band Hubbard Hamiltonian is given in standard notations as 

T-(- = tdpJ2 i4aPja + h.c.) + tpp J2 ip]aPj'^ + h.c.) 

+ AeY.n% + UdY.nf^nf^ (1) 

ja i 

where d^ creates a Cu3dx2_y2 hole and p^ an 02po- hole, tdpitpp) is the nearest-neighbor d-p 
{p-p) transfer. Here the repulsion within the p orbitals and the repulsion between d and p 
orbitals have been neglected for simplicity. 

Great efforts have been made to search for superconductivity in this model [0-0], but 
indications of the off-diagonal long-range order have not been detected so far. There is also 
a variational Monte Carlo study []^, but the justification of the variational wave functions 
remains somewhat open. 

Subsequently reductions to effective Hamiltonians as certain limits of the original three- 
band model have been attempted. In the limit of large level offset {Ae ^ Ud,tdp), the 
system is equivalent to the single-band Hubbard model with the on-site interaction Ud and 
the effective nearest-neighbor hopping tes = tdp/^^- If further put U ^ t^s, the system 
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reduces to the t-J model with excluded double occupancies and J = itlff/U. Thus the t-J 
model is a natural limit of the three-band model in the Mott-Hubbard regime. 

However, the real cuprates lie in the charge-transfer regime. Zhang and Rice have 
proposed that even in this case, the low-lying states of the three-band model may be essen- 
tially represented by the t-J model, at least in the limit of Ud ^ Ae ^ t^p, and provided 
that the spin-triplet d-p molecular orbitals may be neglected. The t-J model thus picks 
up spin-singlet d-p molecular orbitals (or bonding Kondo states), which are envisaged to 
experience a superexchange interaction J while moving around with an effective transfer t 
with double occupancies avoided. The superexchange provides a natural source of an effec- 
tive attraction among the molecular orbitals, and extensive theoretical works have indeed 
indicated that the t-J model superconducts for a certain range of J/t. In one dimension 
(ID) this is shown clearly from a phase diagram having a finite pairing-dominated region 



around J ~ 2i(: 0. In 2D, exact diagonalization results [|T0| indicate that the dx2_y2-wave 
paring correlation function is long-tailed for sufficiently large J ^ t, which is also supported 
from variational Monte Carlo studies [|TTH T3|. 

Now, even if the t-J model can be superconductive, the following fundamental questions 
do remain for the original three-band Hubbard model: 

(i) Does the perturbative picture that maps the three-band model into t-J model in the 
limit of tfip/ Ae,tap/Ud —>■ remain valid for finite, realistic values of parameters ? In real 
materials Ae ~ 2.5tdp [0 is only moderate, where the validity of the perturbation is not at 
all clear. 

(ii) Even if the perturbation is to remain valid through e.g. renormalizations, whether the 
resultant J/t can become large enough to guarantee a high Tc is a nontrivial question. 

(iii) Would there be a qualitative difference between the Mott-Hubbard and charge-transfer 
regimes concerning the appearance of superconductivity via e.g. different effective J/t men- 
tioned in (ii) ? 

All these points evoke another basic question, i.e., does the single-band Unhhsnd model, 
which shares the t-J model as an effective Hamiltonian in the strong-correlation limit, have 
a superconducting phase ? In ID, the conformal field theory indicates that no matter 
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how U /t is increased, the superconducting correlation fails to become dominant, |1T6|JT5 



indicating a behavior distinct from the situation when we let J ~ t in the t-J model. In 
the 2D Hubbard model, quantum Monte Carlo calculations up to f/ = 4t still show no sign 
of the off-diagonal long-range order |l7|. To reconcile this, we have to consider a possibility 
that either the effective J/t is small, oi U = At is already outside the perturbative region. 
If the single-band Hubbard model remains normal for the whole range of parameters, while 
the three-band Hubbard model with finite, realistic values of parameters does superconduct, 
the Mott-Hubbard and charge-transfer regimes may possibly belong to different universality 
classes. 

These problems have remained a long-standing puzzle, which is exactly our motivation 
to revisit the three-band Hubbard model, where we cover a hitherto unfathomed range 
of parameters. If the answer is positive, we will have a stronger ground to consider the 
superconductivity in cuprates to be of electronic origin. 

We employ the quantum Monte Carlo (QMC) method, where our motivation described 
above calls for special emphasis upon the following. 

(i) We consider the range of Ae and Ud extending to the bare width, W , of the most relevant 
(Cu3(i-O2po- anti-bonding) band. We define the case where both Ae and Ud are comparable 
with W to be the strong-correlation regime in the following sense. The relevant energy to 
be compared with W should be the effective repulsion within the d-p Wannier orbital, which 
should be greater than MinjAe, Ud]-, the minimum cost of energy for two holes occupying 
the same Wannier orbital. This is in fact illustrated in the low-lying spectra of finite systems, 
where the levels of the three-band model with Ae: = 3.6eV and Ud = lO.SeV are best-fit to 
those of the single-band Hubbard model with U ~ 5eV [Q. 

(ii) The carrier doping is kept close to the experimentally known optimum value {5 ~ 0.15) 
for the superconductivity. 

(iii) Since a reliable detection of the pairing correlation is required, we adopt the ground- 
state (or projector) QMC formalism with the projection imaginary time of at least 12tdp to 
ensure convergence. 

(iv) The sample-size dependence is studied for lattice sizes up to 8 x 8 unit cells (192 atoms), 
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which is combined with a real-space analysis to probe the range of the pairing correlation. 

To our knowledge, previous calculations do not satisfy all of these conditions simultane- 
ously. For the largest Ae and Ud considered here, the CPU time required was typically 50 
hours on HITAC S-3800 supercomputer. 

As for the symmetry of the pairing, we have considered (ia.2_y2-wave {fd = cosg^ — cos qy) 
and extended s-wave {fs = cosqx + cos qy) pairing, for which we have calculated the k = 
Fourier component of the real space correlation function, Sa = ^^^{Al^Aa + AqAJ^), with 

Aa = Eq/a(q)(c?qTC^-qi +PqTP-qi +PqTP^qi)• 
We first focus on the hole doping. We go from 18 holes for 4x4 unit cells (doping 

6 = 0.125), 42 for 6 X 6 (5 = 0.166), to 74 for 8 x 8 (5 = 0.156). These fillings are chosen 

so as to satisfy (i) the proximity to 5 ~ 0.15, and (ii) the closed-shell condition (with a 

non-degenerate one-electron ground state) to ensure the stability of the QMC calculation. 

We have set tpp = —OAtdp jl^ . 

In Fig.0 the dependence of 5*^ on Ae (a) or Ud (b) is shown. For small Ae and/or Ud, 
Sd decreases with Ae and Ud- An increase in Ae or Ud implies an increased ratio (electron- 
electron repulsion) / (band width) , and the result in the weakly correlated regime shows that 
this works unfavorably for superconductivity as naively expected. However, Sd dramatically 
begins to increase with these parameters for larger values of Ae and/or Ud- The crossover 
to this behavior occurs in the 'strong-correlation' regime where both Ae and Ud exceed the 
band width W of the anti-bonding d-p band (W ~ 2.33tdp for Ae = 2-7tdp and tpp = —O.itdp)- 

Right above the strong-correlation regime the pairing correlation starts to grow with the 
system size, which can be interpreted as a tendency toward the formation of off-diagonal 
long-range order. This is in sharp contrast with the weak-correlation regime, where 5*^ has 
a small, inverse size dependence. 

To check that we are really looking at the long-range part of the pairing correlation, we 
have looked into their behavior in a real space. If we decompose into a sum over the 
real space distance Ar, Sa = Saf ^^(Ar) with Sai^Ar) being the correlation function in real 
space. In Fig.| we represent Srf(Ar) by Sd{R) defined by restricting the sum in the above 
formula to |Aa;|, \Ay\ < R (in the periodic boundary condition) , where Ar = (Ax, Ay). We 
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can see that Sd{R) monotonically increases as we include more and more distant correlations, 
which implies that the growth of the k = component, Sd, is indeed caused by the extension 
of the pairing correlation beyond the system size. 

An indication that this kind of caution is indeed necessary is shown in the inset of Fig.|^. 
Namely, although the extended s-wave pairing correlation, Ss, also increases with the system 
size, its real-space behavior, Ss{R), is almost a constant, indicating that the size dependence 
only signifies a short-range correlation. 

Now, despite the absence of electron-hole symmetry in the three-band Hubbard model, 
the electron-doped materials such as Ndi_2:Ce^Cu04 have been shown to superconduct as 
well. Quite apart from this, the electron-doped case is of another theoretical interest in the 
following sense. In the limit of Ae —>■ oo, Ua oo the electrons are doped to (i.e., holes 
are taken out from) the Cu3d orbital, with only small amount of carriers left in 02p. This 
should make the system closer to the single-band Hubbard model than in the hole-doped 
case. 

In doping electrons (or in taking out holes in the hole picture) the best choice of the band 
filling satisfying the above conditions are 58 holes /8 x 8 (with the doping level 6 = 0.1) and 
26 holes /6 X 6 {6 = 0.26). Thus, the band filhngs are unfortunately not so close for the two 
sizes, which makes the analysis less conclusive. 

Nevertheless, does again become an increasing function of Ae and Ud for Ae ~ 2.5tdp 
and Ud ~ ^-Otdp, where Sd grows as the size becomes 8x8 from 6x6, indicating a tendency 
toward dj;2_y2 electron-pairing in the strongly-correlated regime. 

Encouraged by the electron-doped result, we move on to our final motivation. Namely 
we go back to the hole-doped case to investigate the Mott-Hubbard regime {Ud < Ae) with 
large Ae, which leaves few 02p holes to give another natural way to approach the single-band 
Hubbard model as mentioned earlier. 

In Fig.^ we show the dependence of Sd on Ae with a fixed Ud = l-8tdp (a) or on Ud with 
a fixed Ae = S.Qtdp (b) with the same system sizes and band fillings as in Fig.^ Strikingly 
enough, the system size dependence does appear as well for larger Ae and U just like in 
Fig-E- 
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If we now combine this result with that in the hole-doped charge-transfer regime, the 
following picture emerges. Suppose we compare the relevant energy in the Mott-Hubbard 
regime, MinjAe, Ud} = Ud, with the width of the anti-bonding band, which is W = 1.87tdp 
for Ae = S.Qtdp (and tpp = — O.St^p which we have assumed here). The region at which the 
pairing correlation emerges is precisely Ud ~ W, which is a counterpart to Ae W in the 
charge-transfer regime. Hence, no matter which regime in the three-band model, a tendency 
toward (i^2_j^2-wave pairing superconductivity emerges when the relevant energy {Ud or Ae) 
exceeds W, i.e., when our definition of the strong-correlation criterion is met. In principle 
we can thus envisage a superconductivity phase diagram against Ae and Ud, something like 
UAe > 8t^p if we use the perturbative expression for W. 

If we now recall our reasoning, the result summarized above amounts that either (i) the 
single-band Hubbard model with U/ 1 as large as the bare band width should concomitantly 
exhibit superconductivity, or (ii) we are looking at a regime where the finiteness of Ae makes 
the universality class of the three-band model distinct from that of the single-band Hubbard 
model through e.g. different ranges of the effective J/t. The former possibility that the 
three-band model already resembles the one-band Hubbard model when Ae is increased up 
to ?>Mdp does not contradict with the previous one-band QMC results, where the largest 
U so far studied is only half the band width. At. U = \W is mimicked by the three-band 
model with Ud ~ tpd for Ae = 3.6tdp, for which the sign of pairing is certainly absent in the 
present result as well. We believe this problem deserves further investigations. 

Finally we comment on the possible relevance of our result in the charge-transfer regime 
to the high Tc materials. The value of Ae ~ 2.5tdp where Sd grows with system size is 
remarkably close to the value {Ae = 2.7tdp) obtained from a first-principles calculation 
for La2Cu04 As for the value of Ud, the maximum value tractable with the QMC 

{Ud/tdp = 3 ~ 3.5) happens to be smaller than realistic values {Ud/tdp = 6 ~ 8), but even 
for these small values the tendency for superconductivity already emerges. We expect that 
the tendency can become stronger for larger values of Ud- 

In summary, we have detected an indication of superconductivity in the three-band 
Hubbard model without reducing it into some effective model. Strikingly, this indication 
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emerges in both of the charge-transfer and Mott-Hubbard regimes, and in both of the hole- 
doped and electron-doped cases as long as we are sitting in the strong-correlation regime. 
Since all of these regimes and cases share the t-J model as some limiting cases, this might 
suggest a scenario in which the superconductivity as conceived in the t-J limit extends well 
into the realistic parameter regime. 

Numerical calculations were done on HITAC S3800/280 at the Computer Center of the 
University of Tokyo, and FACOM VPP 500/40 at the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo. This work was supported by 'Project for Parallel 
Processing and Super Computing' at Computer Centre, University of Tokyo, arranged by 
Prof. Y. Kanada, and also by Grant No. 07237209 from the Ministry of Education, Science, 
and Culture, Japan. 
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FIGURES 

FIG. 1. The dr^2_y2-waye pairing correlation, Sd-, is plotted (a) against Ae for a fixed Ud = 3.2, 
and (b) against Ud for a fixed Ae = 2.7. We assume the hopping integrals tdp = 1, tpp = —0.4. 
Number of holes and the sizes of the system are 18 holes/4 x 4 unit cells (A), 42 /6 x 6 (Q)) and 
74 /8 X 8 (□). For 8 x 8 a wider range is displayed in the inset of (a) to show the change in sign 
of the gradient. The dashed lines are guide for the eye. 

FIG. 2. The c?3.2_j^2-wave pairing correlation, Sd{R), and the extended s-wave pairing correla- 
tion, Ss{R) (inset) are plotted as a function of the range, R, in real space. 

FIG. 3. Similar plot for Sd as in Fig.l in the Mott-Hubbard regime, U < Ae. Sd is plotted (a) 
against Ae for a fixed Ud = 1.8, and (b) against Ud for a fixed Ae = 3.6 for tdp = I, tpp = —0.3. 
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